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Theory of carrier density in multigated doped graphene sheets with quantum correction 
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The quantum capacitance model is applied to obtain an exact solution for the space-resolved caiTier density in 
a multigated doped graphene sheet at zero temperature, with quantum correction arising from the finite electron 
capacity of the graphene itself taken into account. The exact solution is demonstrated to be equivalent to the 
self-consistent Poisson-Dirac iteration method by showing an illustrative example, where multiple gates with 
irregular shapes and a nonuniform dopant concentration are considered. The solution therefore provides a 
fast and accurate way to compute spatially varying carrier density, on-site electric potential energy, as well as 
quantum capacitance for bulk graphene, allowing for any kind of gating geometry with any number of gates and 
any types of intrinsic doping. 
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Introduction. Manipulation of carrier density in graphene 
by electrical gating is one of the key techniques for graphene 
electronics. Since the first successful isolation of monolayer 
graphene flakes, conductance (resistance) sweep using a sin- 
gle backgate has been a standard electronic characterization 
tool for graphene.' Double-gated graphene opens possibili- 
ties for experimental investigations of graphene pn and pnp 
junctions,^"^ which allow for exploration of the interesting 
physics of Klein paradox^ in graphene. ^^ In order to improve 
the junction quality, graphene he teroj unctions using contact- 
less topgates^''" and embedded local gates "'^^ were proposed 
and investigated. 

More complicated gating geometry is involved in recent 
proposals for graphene-based devices, such as a switching 
device with two topgates,'^ graphene transistors with self- 
aligned gates made by standard patterning with a regular cross 
section,'^ core-shell nanowires with round cross sections,'^ 
or deposited films with T-shaped cross sections.'^ Transport 
through bilayer graphene with multiple top gates up to eight 
was recently investigated;'^ patterning periodic top gates'^ 
on graphene to form quasi-one-dimensional superlattice is, 
in principle, feasible. Whereas a successful transport simu- 
lation relies decisively on the preciseness of the on-site po- 
tential profile, or equivalently the carrier density profile,'^ a 
more reliable theory to deal with general gating geometry is, 
therefore, imperative. 

The theory of gate-induced carrier density started from 
the simplest classical capacitance model,' which regards the 
graphene-substrate-backgate as a parallel-plate capacitor and 
the relevant carrier density in graphene as the surface charge 
density (divided by electron charge — e) induced by the gate. 
Without taking into account the quantum correction due to 
the finite capacity of graphene itself for electrons to reside, 
this model can be straightforwardly generalized to arbitrary 
gating geometry by treating graphene as a perfect conducting 
plane with fixed zero potential. A more precise computation 
of the gate-induced carrier density, however, needs to take into 
account the relation between the induced charge density on 
graphene and the electric potential energy that those charge 
carriers gain, through the graphene density of states. ^""^^ The 
solution to the carrier density with such a correction taken into 
account requires a self-consistent iteration process "''^^'^'' that 



may be suitably termed the Poisson-Dirac method but actu- 
ally corresponds to the quantum capacitance model,^^ where 
an exact solution for single-gated pristine graphene at zero 
temperature has been derived. ^^ 

In this paper, the spatial profile of carrier density in mono- 
layer graphene due to arbitrary gating and doping is exactly 
solved within the quantum capacitance model. The solution 
has been further tested by comparing with the self-consistent 
Poisson-Dirac method, showing very good agreement be- 
tween the two and, hence, their equivalence. A numerical 
example will be illustrated at the end. Throughout, we will 
restrict our discussion to bulk graphene at zero temperature 
and approximate the energy dispersion within the linear Dirac 
model, E — dnhvyk, which leads to the density of states (per 
unit area) linear in energy, D{E) =2\E\ /n{Ti\'F)- The carrier 
density is given by integrating the density of states over the 
energy. 
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which is the underlying origin of the quantum correction to the 
gate-induced graphene carrier density in the following deriva- 
tions. We are, therefore, working in the single-particle picture, 
and the solution within the quantum capacitance model to be 
presented is exact in the sense that no iteration is required dur- 
ing the solution process, as contrary to the following Poisson- 
Dirac method. 

Self-consistent Poisson-Dirac iteration method. Consider 
a graphene sheet laid in the x-y plain at z = 0. In the presence 
of a dopant concentration nQ{x,y) without electric gating, the 
quasi-Fermi level is given by 



£'o(jc,y) = ?,gn[nQ{x,y)]hvF yj % |no(-«,3') | , 



(2) 



which is obtained from Eq. (1). When gate voltages of, in gen- 
eral, A^ metalic gates are applied as sketched in Fig. 1(a), the 
electron in the graphene layer at (x,^) gains an electrostatic 
potential energy — eVG(x,y), where — e is the electron charge 
and VG(.x:,y) = «(x,y,0) is the electrostatic potential m(x,3',z) 
at z = to be numerically solved from the Poisson equation 
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FIG. 1. (a) Schematic of a graphene sheet subject to A^ metalic 
gates, (b) Equivalent circuit plot of (a) with quantum capacitance 
of graphene Qc taken into account. 



with Eq the permittivity in free space and er{x,y,z) the relative 
permittivity that can be in principle position dependent. 

The energy gain of the electron implies the raising of the 
energy band of graphene and, hence, the lowering of the quasi- 
Fermi level. The graphene carrier density n therefore obeys 
Eq. (l)with 

E(x,y) ^ Eo{x,y) - [-eVG{x,y)] 
Tivp hvf 

= sgn[«o(^,3')]v/^rK(^+^^|^, (4) 

where EQ{x,y) is given by Eq. (2). Together with the charges 
of the dopant ions that maintain the neutrality of the graphene 
sheet, the net charge density on graphene divided by Eq is 
given by 



p{x,y) 



eo 



= - { «o(x,y) --^gn[E{x,y)]E{x,yf \ , (5) 



where E{x,y) = E{x,y)/hvf is given by Eq. (4). Equation (5) 
is the boundary condition at the graphene sheet for the Poisson 
equation (3). This boundary condition contains the solution 
'^G{x,y) =u{x,y,z = Q) and, hence, makes the solution process 
iterative. 

Quantum capacitance model. The system of A^ metalic 
gates labeled by j ~ 1,2,--- ,A^ plus the graphene sheet la- 
beled by G as sketched in Fig. 1(a) is equivalent to the circuit 
plot shown in Fig. 1(b), where the quantum capacitance of the 
graphene sheet Cq is considered. Regarding G as the refer- 
ence conductor with electric potential Vq, the charge density 
on the surface of each gate can be expressed as 

Pi - Cig(Vi - Vg) + Cn{Vx - V2) + • • ■+Cm{Vx - Vn) 
Pi = Cn{V2 - Vi) + C2g(V2 - Vg) + • • • + C2n{V2 - Vn) 



density on G should be the negative of the total charge den- 
sity on the A^ metalic gates: pc = —Y!j=iPj- The net elec- 
tron number density on G is, therefore, hg = pc/i—e) = 
Y!j=iCjGiVj —V(})le. Suppose there is an intrinsic doping 
concentration of «o in graphene. The net charge density on 
G is not affected since the number of doped electrons should 
equal the number of dopant ions, Pc ^> Pg + ^"0 ~ e«o == Pg- 
The net carrier density of graphene, however, is given by 
n = (pq — eno)/{—e) = n^ + nQ, which should obey Eq. (1), 
i.e., nc + no = sgn{EQ + eVc)[{Eo + eVc)/J'iVF]^/7^, just like 
in the Poisson-Dirac method. We therefore need to solve the 
quadratic equation for Vg, 
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(7) 
After some tedious but straightforward algebra, the carrier 
density of graphene in the presence of dopant concentration 

«o and A^ gates with voltages Vi , ■ ■ ■ , Va? is given by 



n = nc + sgn(nc)ne I 1-^/1+2 ) +sgn(«o)-y/2«Q |no|, 

(8) 




where 



N r 



is the classical contribution from doping and gating, and 



(9) 
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arises solely from the quantum capacitance, leading to the 
second and third terms in Eq. (8) as the quantum correction. 
Equations (8)-(10) with A^ = l,«o = 0, and «c > clearly re- 
cover the results for single-gated pristine graphene given in 
Ref. 22. Contrary to the undoped case,^^ the third term in Eq. 
(8) is responsible for the shift of the quasi-Fermi level due to 
doping and is typically weak for a reasonable «o. 

In addition to the doping concentration «o that can have 
any kind of spatial profile, the position dependence enters 
the carrier density (8) through the self-partial capacitances 
CigiQgi" ,Cng, which can be computed numerically but 
exactly. For the /th gate, by grounding all the other conduc- 
tors including the graphene sheet, i.e., Vj^i = and Vg = 0, 
Eq. (6) suggests «c = -I^=i Pj/{-e) = {CiG/e)Vi. The self- 
partial capacitance for gate / is, therefore, given by 



PN = CwiVN -Vi)+ C2NiVN - V2) + • • • + Cng{Vn ~ Vg) 

(6) 

where Cig, ■ ■ ■ ,Cng are self-partial capacitances and Cij with 
i 7^ j are mutual partial capacitances.^^ Since the whole iso- 
lated system should remain charge neutral, the net charge 
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where «c = ±£r£o(<5M/<5z),=o±/e can be numerically com- 
puted by any kind of finite-element simulator 

With the definitions (9) and (10), one may also write the 
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solution Vg to Eq. (7), 



Vg 



sgn(«c)nG(l-Wl+2^) +sgn 



(n())V2«e|no 
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(12) 

which has a reasonable form of charge divided by capacitance, 
with the numerator containing only the quantum coiTection 
terms in Eq. (8). The absence of nc in the numerator of Vg 
agrees with our earlier remark that the classical capacitance 
model regards graphene as a perfect conducting plane with 
fixed zero potential so «c does not contribute to Vq. 

Equation (12) allows for a direct comparison with the iter- 
ative solution obtained from the self-consistent Poisson-Dirac 
method, as we will show with an explicit example soon. Mul- 
tiplying Eq. (12) with the electron charge together with the 
quasi-Fermi level shift £0 due to doping, —(£0 + eVo) pro- 
vides for the graphene transport calculation a realistic on-site 
energy profile that guarantees a reliable quantum transport 
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FIG. 2. (Color online) (a) Side view of a graphene sheet (with a 
hyperbolic-tangent-shaped intrinsic doping iiq) and a backgate (with 
Kbg = —20 V) sandwiching a SiOa with two embedded local gates 
(with Vigi = —1.8 V and Vig2 = 1.5 V); the color shading shows 
the electric potential u{x,z) obtained by the self-consistent Poisson- 
Dirac method, (b) The electric potential at the graphene layer V(j{x) 
obtained by the Poisson-Dirac method and the quantum capacitance 
model. 



simulation; see, for example, Ref. 19 for the case with ne- 
glected «o. Furthermore, the channel electrostatic potential 
Vg given in Eq. (12) also allows us to write down the quan- 
tum capacitance of the graphene sheet in the low-temperature 
limit:22 Cq « (2/7r) {e/hvpf leVcl- 

Numerical example. Armed with the above introduced 
theories, we next numerically demonstrate the equivalence of 
the quantum capacitance model to the self-consistent Poisson- 
Dirac iteration method by considering a specific example. To 
be simple but general, let us consider a quasi-one-dimensional 
system along x with translation invariance along y, composed 
of a doped graphene sheet gated by one flat backgate and two 
embedded local gates with irregular shapes roughly 10 nm un- 
der graphene; see Fig. 2(a). Embedding such local gates at 
such a shallow depth allows independent control of the car- 
rier density in the locally gated region due to screening of 
the backgate contribution and can be experimentally achieved; 
see, for example, Ref. 1 1 . The finite-element method is imple- 
mented in the iteration process for the Poisson-Dirac method 
as well as the exactly solvable self-partial capacitances [Eq. 
(11)] for the quantum capacitance model, and the pdetool in 
MatlaB'^ is chosen as the simulator for the present demon- 
stration. 

The electric potential ii{x,z) shown in Fig. 2(a) is obtained 
by the self-consistent Poisson-Dirac method with backgate 
voltage Vbg = —20V and local gate voltages Vigi = —1.8V 
and V\a2 = 1 .5 V and an intrinsic doping described by no(x) = 
—5 X 10" tanh(.ic/40)cm^2, where the position coordinate x 
is in units of nm. The iterated potential solution Vg{x) = 
u{x,z = 0) at the graphene layer is compared in Fig. 2(b) with 
the exact solution (12) obtained within the quantum capaci- 
tance model, showing an excellent agreement with each other. 
With other gate voltages and other shapes of no{x), the agree- 
ment remains exact. Note that the numerical example chosen 
here is basically a complicated version of Ref. 11, including 
the proper range of the gate voltages, except that an artificial 
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FIG. 3. (Color online) Carrier density profiles of intrinsic doping 
n()(x), classical capacitance model ncc(^)> Poisson-Dirac method 
«pd(x), quantum capacitance model iiqc(x), and the difference 
npD(x) — «Qc(^), with identical parameters used in Fig. 2. 
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doping profile «o with hyperbolic tangent shape is considered, 
in order for the comparison to be general. 

The spatial profiles of the carrier densities n(j{x), ncc{x), 
npoC'c), iiQcix), as well as the difference np£>{x) —nQc{x) are 
shown in Fig. 3. Here the subscripts CC, PD, and QC de- 
note "classical capacitance," "Poisson-Dirac," and "quantum 
capacitance," respectively. The carrier density within the clas- 
sical capacitance model «cc is obtained by first computing 
the induced surface charge at z = 0^ with the graphene layer 
grounded (Vg = 0) and then adding the dopant concentration 
«o or, equivalently, by Eq. (9) with the self-partial capaci- 
tances [Eq. (11)] numerically computed. 

As the quantum correction, i.e., the second and third terms 
in Eq. (8), always reduces the magnitude of the net contribu- 
tion of the gates, the classical solution always overestimates 
the gate-induced carrier density. This correction is especially 
salient when the gate is close to the graphene sheet, as is 
clearly observed by comparing np£>{x) or nQc{x) with ncc{x) 
in Fig. 3. In addition, the surface roughness of the embedded 
local gates considered here with such a short distance to the 
graphene sheet (roughly lOnm) further introduces a strongly 
fluctuating potential profile [Fig. 2(b)] as well as the corre- 
sponding carrier density profile (Fig. 3) at the locally gated 
regions. 

As in the case of Vg{x) compared in Fig. 2(b), the agree- 
ment between np£,{x) and «qc(-*^) is rather satisfactory. In Fig. 



3, the discrepancy between the Poisson-Dirac method and the 
quantum capacitance model becomes relatively obvious near 
positions where the surface charge density of the boundary 
condition (5) is changing its sign. This implies that the dis- 
crepancy may stem from the inherent numerical limitation of 
the chosen nonlinear partial differential equation solver. 

Conclusion. In conclusion, an exact solution for the 
space-resolved carrier density in multigated doped graphene 
sheets within the quantum capacitance model has been de- 
rived. With an illustrative quasi-one-dimensional example, 
the exact solution is shown to be equivalent to the self- 
consistent Poisson-Dirac iteration method. The solution there- 
fore provides a fast and accurate way to compute spatially 
varying carrier density, on-site potential energy (key input for 
quantum transport simulation), as well as quantum capaci- 
tance for bulk graphene, allowing for any kind of gating geom- 
etry and any types of intrinsic doping. Moreover, the contact 
doping'^'^^ and its corresponding screening potential^" can 
as well be treated by the presented solution, which therefore 
takes care of all three types of doping in graphene — electric, 
chemical, and contact-induced — in a unified manner 
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